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We present spatially resolved measurements characterizing the stagnation layer between two 
obliquely merging supersonic plasma jets. Intra-jet collisionality is very high (Xu <C 1 mm), but 
the inter-jet ion-ion mean free paths are on the same order as the stagnation layer thickness (a few 
cm). Fast-framing camera images show a double-peaked emission profile transverse to the stagna- 
tion layer, with the central emission dip consistent with a density dip observed in the interferometer 
data. We demonstrate that our observations are consistent with collisional oblique shocks. 



Colliding plasmas have been studied in a variety of con- 
texts, e.g., counterstreaming laser-produced plasmas sup- 
porting hohlraum design for indirect-drive inertial con- 
finement fusion [H-El, forming and studying astrophysi- 
cally relevant shocks [4-8], and for applications such as 
pulsed laser deposition |9[ and laser induced breakdown 
spectroscopy [10]. Physics issues arising in these stud- 
ies include plasma interpenetration [Ill4l6| , shock forma- 
tion [17], and the formation and dynamics of a stagna- 
tion layer between the colliding plasmas [l8l-[20|. In this 
work, we present experimental results of two obliquely 
merging supersonic plasma jets. These and other re- 
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cent jet merging experiments 
to explore the feasibility of forming imploding spherical 
plasma liners via an array of merging plasma jets [23l-[27|. 
which could have applications in forming cm-, /is-, and 
Mbar-scale plasmas for fundamental high energy density 
physics studies 12 81 and as a standoff driver for magne- 
toinertial fusion [23. l29l - [32| . Prior experiments studying 
the stagnation layer between colliding laser-produced or 
wire-array z-pinch [33| plasmas were on smaller spatial 
scales (mm or smaller) that could not be fully resolved 
by measurements. New results in the present work are 
the experimental identification and characterization of a 
few-cm thick stagnation layer between colliding plasmas, 
and the demonstration that our observations are consis- 
tent with hydrodynamic oblique shock theory [341 ]. 

Experiments reported here are conducted on the 
Plasma Liner Experiment (PLX) [35], in which two su- 
personic argon plasma jets are formed and launched 
by plasma railguns [36| . Plasma jet parameters at the 
exit of the railgun nozzle (peak n e ~ 2 x 10 16 cm -3 , 
peak T e « 1.4 eV, Vj et ~ 30 km/s, Mach number 
M = Vj e t /C s jet ~ 14, diameter = 5 cm, and length 
« 20 cm) and their evolution during subsequent jet prop- 
agation have been characterized in detail [35]. The jet 
magnetic field inside the railgun is ~3 T, but the classi- 
cal magnetic diffusion time is a few /is [35[, an d thus 
we ignore the effects of a magnetic field by the time 
of jet merging (> 20 /is). Experimental data are from 
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FIG. 1. Schematic of the experiment showing the spheri- 
cal vacuum chamber, location of railgun nozzles mounted 24° 
apart, two obliquely merging plasma jets, (R,Z) coordinates 
used in the paper, approximate interferometer/spectrometer 
lines-of-sight [Z ~ 84 cm), and CCD camera field-of-view. 



an eight-chord laser (561 nm) interferometer [37|, |38|, 
a visible-to-near-infrared 0.275 m survey spectrometer 
(600 lines/mm grating and 0.45 /is gating on the 1024- 
pixel MCP array detector), and an intensified charge- 
coupled device (CCD) visible-imaging camera (DiCam 
Pro, 1280 x 1024 pixels, 12-bit dynamic range). The in- 
terferometer and spectrometer chords intersect the merg- 
ing jets at Z w 84 cm (Fig. [I]). Each interferometer probe 
beam has a diameter of ~ 3 mm, and the chords are each 
separated by 1.5 cm transverse to the merge plane. The 
spectrometer has an w 7 cm diameter field-of-view in the 
vicinity of the merge plane. More details about PLX, the 
railguns, diagnostics, and characterization of individual 
plasma jets are given elsewhere [35|. 

Figure [2] shows the time evolution of oblique jet- 
merging. Formation of a stagnation layer between the 
two jets and its double-peaked emission profile in the 
transverse (R) direction are clearly visible. We made 
interferometer and spectrometer measurements of the 
stagnation layer at Z « 84 cm (Fig. [3^) for the cases 
of top-only, bottom-only, and both jets firing. Fig- 
ure [3b shows the interferometer phase shift Ac/) versus 
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FIG. 2. False-color, slightly cropped, CCD images (log inten- 
sity, 20 ns exposure) of oblique jet-merging (shots 1130, 1128, 
1125, 1120, 1134, 1138). The two railgun nozzles (« 46 cm 
apart) are visible on the right-hand-side of each image. 
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time at the R = 2.25 cm chord position for each case. 
Merged-jet measurements show that, at R = 2.25 cm, 
Emerge > A(j) t0 p + A(j) boUo m (Fig. [3b), implying that 
simple jet interpenetration cannot account for the ob- 
served A(j) of the merged-jet stagnation layer (more quan- 
titative analysis given later). However, at large R (e.g., 
6.75 cm), A^merge & A(j) t0 p + ^bottom (not shown), con- 
sistent with jet interpenetration. Figure [3b shows Acj) 
vs. R at four times for a merged jet. The A<j) dip at 
R = 0.75 cm and peak at R = 2.25-4 cm are well- aligned 
with the emission dip and peak (Fig. [3^), respectively. 

The interferometer A<j) measurements are used to es- 
timate the ion plus neutral density n to t- Our interfer- 
ometer is sensitive to bound and free electrons in the 
plasma [IH, so Acj) contributions from all species and 
ionization states must be considered. The A(j) satisfies 
J ritotdl = c [z^-Err] > wnere the integral is over the 
chord length, C e = Ae 2 /47reom e c 2 is the phase sensi- 
tivity to electrons, Z e ^ = n e /n to t is the mean charge, 
Err = k (2ir / \C e n to t)Kj,k m k n j,k is the error in the 
phase shift due to all ionization states j of all ion species 
/c, rrik is the atomic mass of ion species /c, and Kj^ is the 
specific refractivity. Uncertainty in plasma jet composi- 
tion (due to impurities) can be accounted for by bound- 
ing rttot at the extremes Err = and Err = Err max = 
(27rra/ ' \C e )K rnaxi where K max is the largest specific re- 
fractivity of all the species present. Since we do not 
know precisely the impurity fraction or mixture ratios 
of impurities, we perform our data analysis by consid- 
ering two extreme cases of (i) 100% argon and (ii) 30% 
argon with 70% impurities. The latter is chosen based 
on the difference in measured chamber pressure rise for 
gas-injection-only versus full plasma discharges. Iden- 



FIG. 3. (a) CCD line-out, interferometer chord positions, and 
spectrometer view overlayed on a CCD image (t = 38 /xs). 
(b) Multi-shot averaged interferometer phase shift vs. time 
at R = 2.25 cm. (c) Phase shift vs. chord position for a 
jet-merging case (shot 1120). 



tification of bright Al and O spectral lines in our data 
suggest that impurities are from the alumina-based rail- 
gun insulators. Thus, for case (ii), we approximate the 
plasma jet to be 43% O and 24% Al and assume that 
Err max = Err ai i- Top-jet-only experiments provide 
single-jet Acj) vs. time at Z w 84 cm and R = 2.25 cm 
(Fig. [3b). The average single-jet peak phase shift is 
A<j) = 4.3 ± 0.3° for the data considered (shots 1265- 
1267). All chord positions R = 0.75-11.25 cm (Fig. Hi) 
have similar A</> w 4°. Using Acj) = 4°, Z eff = 0.94 
(inferred from spectroscopy analysis [35[ assuming 100% 
argon), Err max = 0.082, and a jet diameter of 22 cm 
(from CCD images [35|), gives a single-jet density range 
of ntot = Ttsingie = 2.1-2.3 x 10 14 cm -3 . The latter result 
changes by only a few percent for the 30%/70% mixture, 
which changes the inferred Z e R to 0.92. 

By using the interferometry estimates for merged-jet 
density (shown later) and comparing our spectroscopy 
data with non-LTE spectral calculations using Prism- 
SPECT [39[ , we infer Z e ff and T e of the stagnation layer 
at Z w 84 cm. PrismSPECT results for Z eff and T e are 
sensitive to the specific plasma mixture used. Based on 
the presence of certain Ar II lines in the data and by 
comparing to PrismSPECT results, we bound estimates 
of Z eff and T e using the 100% argon and 30%/70% mix- 
ture cases. For the former (Fig. |4^), we infer that peak 
T e > 1.4 eV and Z eff = 0.94. For the latter (Fig.|4]o), we 
infer that 2.2 eV< peak T e < 2.3 eV and Z eff = 1.3-1.4, 



3 



(a) 



1 

? 0.8 

ro.6- 



Shot # 1122 

-n=8x10 14 cm~ 3 , T =1.4 eV 

e 

PrismSPECT non-LTE 



Ar-ll that indicates T >1 .4 eV 



Z = 90 cm 



450 460 470 480 490 500 510 
Wavelength [nm] 



1 

,.0.8 



(b) 



Shot #1122 
— n=5x10 14 cm -3 , T=2.2 eV 

PrismSPECT non-LTE 
0.3 Ar + 0.43 O + 0.24 Al 



Ar-ll that indicates T >2.2 eV 



450 460 470 480 490 500 510 

Wavelength [nm] 



FIG. 4. Spectral data (gray) and non-LTE PrismSPECT cal- 
culations (black) for the merged-jet stagnation layer (Z « 
84 cm, t = 36 jjls). The PrismSPECT calculations are for 
(a) 100% argon and (b) 30%/70% mixture. Lower bounds on 
peak T e are inferred based on the presence of the Ar n lines 
indicated by asterisks. 



with the upper bounds determined by the absence of an 
Al ill line in the data. For the mixture, M « 9 compared 
to M ~ 14 for pure argon. 

With estimates of the stagnation layer Z e s in hand, 
we estimate the stagnation layer density and compare 
it with the single-jet (un-shocked) density. At R = 
2.25 cm, the average peak = 14.3 ± 2.4° (Fig. [3]o) 
(shots 1117-1196). Using A</> = 14°, chord path length 
of 22 cm, and Z e R — 0.94 (100% argon case), n to t — 
Emerged = 7.5-8.2 x 10 14 cm -3 . In this case the den- 
sity increase n mer g ed /n sin gi e = 3.2-3.8. For Z eff = 1.4 
(30%/70% mixture case), the stagnation layer density is 
ntot = ^merged — 5.0-5.3 x 10 14 cm -3 , and the density 
increase is rtm erge dln S i n gi e = 2.1-2.5. Thus, the observed 
range of n merge d/n S i n gi e = 2.1-3.8, exceeding the factor 
of two expected for jet interpenetration. 

We compare the experimentally inferred density jumps 
with oblique shock theory At Z « 84 cm and 

M = 9-14, the theory predicts a shock angle f3 « 19°- 
20°, as discussed later. For 7 = 1.4, the predicted density 
jump across an oblique shock [28] is n shock /n unshocked = 
(Msin/3) 2 (7+l)/[(Msin/3) 2 (7-l)-2] « 4.0-4.9. Differ- 
ence between the measured and predicted density jumps 
could be due to 3D (e.g., pressure-relief in the out-of- 
page dimension) and/or equation-of-state (e.g., ioniza- 
tion [22j| ) effects not modeled by the theory. 

The stagnation layer thicknesses as observed in the 
merged-jet emission (Fig. [5} and Acj) vs. R profiles 
(Fig. [3t) are similar in scale (few cm). In a collisional 
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FIG. 5. CCD image line-outs versus R (transverse to stagna- 
tion layer) at Z = 90 cm (horizontal pixel # 654 out of 1024), 
corresponding to the images of Fig. [2j Progressive times are 
shown with increasing count offsets to avoid trace overlap. 



plasma, the layer thickness is expected [14] to be of or- 
der the counter-streaming ion-ion mean free path (mfp) 
Xu> ~ v re i/4vu> [22[, where v re i is the relative transverse 
velocity between obliquely merging jets and the slowing- 
down rate in the fast approximation (v re i ^> vu) is 
v iV « 9.0 x 10- 8 (l/M+l/M')(M 1/2 A 3/2 K'(^eff^ ff ) 2 In A 
[40[. Note that in our parameter regime, the inter-jet 
\e ^ Xu- We estimate Xu> by considering jets of 100% 
argon and the 30%/70% mixture previously discussed, 
in all cases using v re \ = 20 km/s. For Ar-Ar stopping, 
Xu = 3.47 cm (for m = 8 x 10 14 cm" 3 , T e = 1.4 eV, 
Z eff = 0.94). Pure Al-Al and 0-0 stopping yield Xu = 
0.16 and 0.62 cm, respectively (for rti = 5 x 10 14 cm -3 , 
T e = 2.2 eV, Z e ff 5 Ai = 2.0, Z e ^^o = 1.0). For inter-species 
collisions in a mixed-species jet, using the mixture given 
in Fig. [3)3, Xu/ w 0.57-6.18 cm. We have also estimated 
the inter-jet mfp due to Ar + -Ar charge-exchange and mo- 
mentum transfer [4l| to be ~ 3 cm. These estimates im- 
ply that our inter-jet merging is in a collisional regime, 
which is consistent with a more detailed treatment of 
inter-jet ion- ion stopping including jet profile effects (22| . 

Assuming that the emission layers are post-shocked 
plasma, we postulate that their edges (at larger \R\) are 
the shock boundaries (Fig. [6]). Qualitatively, the merg- 
ing geometry for an individual jet resembles that of a su- 
personic flow past a wedge [34[ and should result in the 
formation of a linear attached shock at an angle j3 with re- 
spect to the original flow direction. The angle /? is a func- 
tion of M and the wedge angle 8 between the jet flow and 

M^sin 2 ^-! i 



the midplane, and satisfies tan 6 = 2 cot j3[ 



M 2 (7+cos2/3)+2J 



[34|. In our case, tanc) = (23 cm)/Z^, where Z{ is the 
point at which jets first interact. We estimated Zi from 
CCD images as the minimum Z for which merged-jet 
emission is observed, with Zi « 45 cm at t = 26 /is 
evolving to Zi w 25 cm at t = 46 /is. For this range 
of Zi, we predict shock angles (relative to the midplane) 
P - S « 8°-17° (M = 9) and 7°-15° (M = 14). For 
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FIG. 6. CCD image (shot 1089, t = 30 fis) with postulated 
shock boundaries (solid white lines) and observed shock angle 
/3 — 5 ~ 5° relative to the midplane. 



shot 1089 with Z { « 30 cm, f3 - S w 5° (Fig. [6]), which 
is within a factor of two of the predicted j3 — 5 ~ 9°- 
10°. This is reasonable given that the prediction does 
not include 3D nor equation-of- state [22| effects. Oblique 
shock theory also predicts the formation of a detached 
shock [34[, which would have a curved shock boundary 
for S > S max « 45° (for M = 9-14 and Z { « 25 cm 
at late times). CCD images at t > 42 /is (Fig. [2j) show 
curvature of the emission layers away from the midplane, 
suggestive of detached shocks. 

To further evaluate the consistency of the experimen- 
tal results with oblique shock theory [34[, we ran ID 
multi-fluid simulations of merging jets, treating the elec- 
trons as one fluid and the ions of each jet as a second 
and third fluid. This models the transverse (R) dy- 
namics of the oblique merging. Simulations were per- 
formed with the multi-fluid magnetohydrodynamics code 
USim, formerly known as Nautilus [42|. In the simula- 



tions, the jets are assumed to be 100% Ar II with ini- 
tial n e = rii = 10 14 cm" 3 , T e = Ti = 1.4 eV, and 
velocities of ±6.2 km/s (i.e., transverse component of 
Vj e t ~ 30 km/s). We used a density profile in the lead- 
ing edge of the jet as shown in the top panels of Fig. 
We performed both collisional [3] and collisionless sim- 
ulations with the same initial conditions to evaluate the 
differences between the two cases. The simulations have 
resolution ~ 100 /im, which does not resolve the Debye 
length (Ad ~ 1.0 /im). However, we verified the for- 
mation of an electrostatic collisionless shock on reduced 
density simulations that do resolve Xp. Details of the 
semi-implicit numerical algorithm are described in [43|. 
At 1 /is after merging begins, both collisionless and col- 
lisional simulations show a small initial density buildup 
at the merging interface. Within 5 /is, a density dip ap- 
pears at the midplane, similar to the structure seen for 
tens of microseconds in the experiments. At these jet 
parameters, the collisional and collisionless simulations 
give the same qualitative features, i.e., the multi-layer re- 
flected shock structure that is consistent with the double- 
peaked emission seen in Fig. [2l Subtle differences be- 
tween the collisional and collisionless simulation results, 
i.e., reflected shock speed and degree of ion interpenetra- 
tion, cannot at present be resolved by the experimental 
data. The simulation results provide further support to 
the interpretation that our observations are consistent 
with collisional oblique shocks. 
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FIG. 7. Density versus spatial dimension from ID multi- 
fluid collisional (left column) and collisionless (right column) 
plasma simulations that model the transverse (R) dynamics of 
our experiments, illustrating the formation of reflected shocks 
in both cases. 



In summary, we have experimentally characterized the 
stagnation layer between two obliquely merging super- 
sonic plasma jets. The jets are individually very highly 
collisional, but the inter-jet ion-ion collisional mean free 
path is of the same order as the stagnation layer thickness 
of a few centimeters. CCD images show the formation of 
a stagnation layer with a double-peaked emission profile 
transverse to the layer, with the central emission dip con- 
sistent with a density dip observed in the interferometer 
data. The overall geometry of the observed stagnation 
layer structure is consistent with hydrodynamic oblique 
shock theory. Furthermore, collisional ID multi- fluid 
plasma simulations that model the transverse dynamics 
of the oblique merging show the formation and evolution 
of reflected shocks with a central density dip consistent 
with the observed stagnation layer. Ongoing experiments 
are now employing lower jet density and higher jet veloc- 
ity to study head-on jet merging with very low inter-jet 
ion-ion collisionality. 
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